
		**********************************************************
		** File: EBA pdest 16-30.do								**
		** Paper: Human Trafficking Indicators: A New Dataset	**
		** Author: Richard Frank								**
		** Date: July 11, 2021									**
		** Task: Running 2nd series of pdest EBA models			**
 		**********************************************************
				clear
				version 16.1
				set max_memory 32g, permanently
				set seed 1234
				set level 95 		
				  cd "C:\Users\Administrator\Desktop\"
				use "HT EBA merged DV and IVs.dta"
	
*******************	
 ** 16. homicides **
 ***************
		
		probit  F.pdest  /**/ homicides /**/
		regsave using /**/ eba_pdest_16 /**/, replace
		
		tuples 																///
		  agr america_l america_n conflict_ext conflict_int 				///
		 cpi easia_pacific energyuse ethnic_tension ethnic  				///
		 europe execleft execright fdi food      							///
		/**/  infant KOFSoGI language legal_british 						///
		 legal_french legal_german legal_scan literacy 						///
		 lp_catho80 lp_protmg80 media mena migration 						///
		 mortality_under5 lp_muslim80 oecd peacekeepersfl peacekeepersl 				///
		 polity2 pop prevent prostitutionlaw protect 						///
		 pve refugees religion rle sasia 									///
		 ssa tourism trade unemployment urban 								///
		 vae visa wecon  english 										///
		 french german spanish portuguese landlocked 						///	
				, max(3) min(3)  kronecker
		
		forvalues i = 1/`ntuples'   {
		 probit  F.pdest  /**/ homicides /**/ `tuple`i'' gdppcl, robust cluster(ccode) 
		regsave using /**/ eba_pdest_16 /**/, tstat pval ci autoid append detail(all) 	
						}							
												
 ***************
 ** 17. infant **
 ***************
		
		probit  F.pdest  /**/ infant /**/
		regsave using /**/ eba_pdest_17 /**/, replace
		
		tuples 																///
		  agr america_l america_n conflict_ext conflict_int 				///
		 cpi easia_pacific energyuse ethnic_tension ethnic  				///
		 europe execleft execright fdi food      							///
		 homicides /**/   KOFSoGI language legal_british 					///
		 legal_french legal_german legal_scan literacy 						///
		 lp_catho80 lp_protmg80 media mena migration 						///
		 mortality_under5 lp_muslim80 oecd peacekeepersfl peacekeepersl 				///
		 polity2 pop prevent prostitutionlaw protect 						///
		 pve refugees religion rle sasia 									///
		 ssa tourism trade unemployment urban 								///
		 vae visa wecon  english 										///
		 french german spanish portuguese landlocked 						///	
				, max(3) min(3)  kronecker
		
		forvalues i = 1/`ntuples'   {
		 probit  F.pdest  /**/ infant /**/ `tuple`i'' gdppcl, robust cluster(ccode) 
		regsave using /**/ eba_pdest_17 /**/, tstat pval ci autoid append detail(all) 	
						}							
												
												
 ***************
 ** 18. KOFSoGI **
 ***************
		
		probit  F.pdest  /**/ KOFSoGI /**/
		regsave using /**/ eba_pdest_18 /**/, replace
		
		tuples 																///
		  agr america_l america_n conflict_ext conflict_int 				///
		 cpi easia_pacific energyuse ethnic_tension ethnic  				///
		 europe execleft execright fdi food      							///
		 homicides infant /**/    language legal_british 					///
		 legal_french legal_german legal_scan literacy 						///
		 lp_catho80 lp_protmg80 media mena migration 						///
		 mortality_under5 lp_muslim80 oecd peacekeepersfl peacekeepersl 				///
		 polity2 pop prevent prostitutionlaw protect 						///
		 pve refugees religion rle sasia 									///
		 ssa tourism trade unemployment urban 								///
		 vae visa wecon  english 										///
		 french german spanish portuguese landlocked 						///	
				, max(3) min(3)  kronecker
		
		forvalues i = 1/`ntuples'   {
		 probit  F.pdest  /**/ KOFSoGI /**/ `tuple`i'' gdppcl, robust cluster(ccode) 
		regsave using /**/ eba_pdest_18 /**/, tstat pval ci autoid append detail(all) 	
						}							 					
									
 ***************
 ** 19. language **
 ***************
		
		probit  F.pdest  /**/ language /**/
		regsave using /**/ eba_pdest_19 /**/, replace
		
		tuples 																///
		  agr america_l america_n conflict_ext conflict_int 				///
		 cpi easia_pacific energyuse ethnic_tension ethnic  				///
		 europe execleft execright fdi food      							///
		 homicides infant KOFSoGI /**/     legal_british 					///
		 legal_french legal_german legal_scan literacy 						///
		 lp_catho80 lp_protmg80 media mena migration 						///
		 mortality_under5 lp_muslim80 oecd peacekeepersfl peacekeepersl 				///
		 polity2 pop prevent prostitutionlaw protect 						///
		 pve refugees religion rle sasia 									///
		 ssa tourism trade unemployment urban 								///
		 vae visa wecon  english 										///
		 french german spanish portuguese landlocked 						///	
				, max(3) min(3)  kronecker
		
		forvalues i = 1/`ntuples'   {
		 probit  F.pdest  /**/ language /**/ `tuple`i'' gdppcl, robust cluster(ccode) 
		regsave using /**/ eba_pdest_19 /**/, tstat pval ci autoid append detail(all) 	
						}							 					
									
 ***************
 ** 20. legal_british **
 ***************
		
		probit  F.pdest  /**/ legal_british /**/
		regsave using /**/ eba_pdest_20 /**/, replace
		
		tuples 																///
		  agr america_l america_n conflict_ext conflict_int 				/// 1-5
		 cpi easia_pacific energyuse ethnic_tension ethnic  				/// 6-10
		 europe execleft execright fdi food      							/// 11-15
		 homicides infant KOFSoGI language /**/         					/// 16-19
		 legal_french legal_german legal_scan literacy 						/// 20-23
		 lp_catho80 lp_protmg80 media mena migration 						/// 24-28
		 mortality_under5 lp_muslim80 oecd peacekeepersfl peacekeepersl 	/// 29-33
		 polity2 pop prevent prostitutionlaw protect 						/// 34-38
		 pve refugees religion rle sasia 									/// 39-43
		 ssa tourism trade unemployment urban 								/// 44-48
		 vae visa wecon  english 											///
		 french german spanish portuguese landlocked 						///	
				, max(3) min(3)  kronecker cond(!(20&27&28))
		
		forvalues i = 1/`ntuples'   {
		 probit  F.pdest  /**/ legal_british /**/ `tuple`i'' gdppcl, robust cluster(ccode) 
		regsave using /**/ eba_pdest_20 /**/, tstat pval ci autoid append detail(all) 	
						}							 					
									
																															
 ***************
 ** 21. legal_french **
 ***************
		
		probit  F.pdest  /**/ legal_french /**/
		regsave using /**/ eba_pdest_21 /**/, replace
		
		tuples 																///
		  agr america_l america_n conflict_ext conflict_int 				/// 1-5
		 cpi easia_pacific energyuse ethnic_tension ethnic  				/// 6-10
		 europe execleft execright fdi food      							/// 11-15
		 homicides infant KOFSoGI language legal_british    				/// 16-20
		 /**/  legal_german legal_scan literacy 							/// 21-23
		 lp_catho80 lp_protmg80 media mena migration 						/// 24-28
		 mortality_under5 lp_muslim80 oecd peacekeepersfl peacekeepersl 	/// 29-33
		 polity2 pop prevent prostitutionlaw protect 						/// 34-38
		 pve refugees religion rle sasia 									/// 39-43
		 ssa tourism trade unemployment urban 								/// 44-48
		 vae visa wecon  english 											/// 49-52
		 french german spanish portuguese landlocked 						///	
				, max(3) min(3)  kronecker cond(!(20&27&28))
		
		forvalues i = 1/`ntuples'   {
		 probit  F.pdest  /**/ legal_french /**/ `tuple`i'' gdppcl, robust cluster(ccode) 
		regsave using /**/ eba_pdest_21 /**/, tstat pval ci autoid append detail(all) 	
						}							 					
									
																																																			
 ***************
 ** 22. legal_german **
 ***************
		
		probit  F.pdest  /**/ legal_german  /**/
		regsave using /**/ eba_pdest_22 /**/, replace
		
		tuples 																///
		  agr america_l america_n conflict_ext conflict_int 				///
		 cpi easia_pacific energyuse ethnic_tension ethnic  				///
		 europe execleft execright fdi food      							///
		 homicides infant KOFSoGI language legal_british    				///
		 legal_french /**/   legal_scan literacy 							///
		 lp_catho80 lp_protmg80 media mena migration 						///
		 mortality_under5 lp_muslim80 oecd peacekeepersfl peacekeepersl 				///
		 polity2 pop prevent prostitutionlaw protect 						///
		 pve refugees religion rle sasia 									///
		 ssa tourism trade unemployment urban 								///
		 vae visa wecon  english 										///
		 french german spanish portuguese landlocked 						///	
				, max(3) min(3)  kronecker
		
		forvalues i = 1/`ntuples'   {
		 probit  F.pdest  /**/ legal_german /**/ `tuple`i'' gdppcl, robust cluster(ccode) 
		regsave using /**/ eba_pdest_22 /**/, tstat pval ci autoid append detail(all) 	
						}							 					
																																			
 ***************
 ** 23. legal_scan **
 ***************
		
		probit  F.pdest  /**/ legal_scan  /**/
		regsave using /**/ eba_pdest_23 /**/, replace
		
		tuples 																///
		  agr america_l america_n conflict_ext conflict_int 				///
		 cpi easia_pacific energyuse ethnic_tension ethnic  				///
		 europe execleft execright fdi food      							///
		 homicides infant KOFSoGI language legal_british    				///
		 legal_french legal_german /**/    literacy 						///
		 lp_catho80 lp_protmg80 media mena migration 						///
		 mortality_under5 lp_muslim80 oecd peacekeepersfl peacekeepersl 				///
		 polity2 pop prevent prostitutionlaw protect 						///
		 pve refugees religion rle sasia 									///
		 ssa tourism trade unemployment urban 								///
		 vae visa wecon  english 										///
		 french german spanish portuguese landlocked 						///	
				, max(3) min(3)  kronecker
		
		forvalues i = 1/`ntuples'   {
		 probit  F.pdest  /**/ legal_scan /**/ `tuple`i'' gdppcl, robust cluster(ccode) 
		regsave using /**/ eba_pdest_23 /**/, tstat pval ci autoid append detail(all) 	
						}							 					
																																			
 ***************
 ** 24. literacy **
 ***************
		
		probit  F.pdest  /**/ literacy  /**/
		regsave using /**/ eba_pdest_24 /**/, replace
		
		tuples 																///
		  agr america_l america_n conflict_ext conflict_int 				///
		 cpi easia_pacific energyuse ethnic_tension ethnic  				///
		 europe execleft execright fdi food      							///
		 homicides infant KOFSoGI language legal_british    				///
		 legal_french legal_german legal_scan /**/     						///
		 lp_catho80 lp_protmg80 media mena migration 						///
		 mortality_under5 lp_muslim80 oecd peacekeepersfl peacekeepersl 				///
		 polity2 pop prevent prostitutionlaw protect 						///
		 pve refugees religion rle sasia 									///
		 ssa tourism trade unemployment urban 								///
		 vae visa wecon  english 										///
		 french german spanish portuguese landlocked 						///	
				, max(3) min(3)  kronecker
		
		forvalues i = 1/`ntuples'   {
		 probit  F.pdest  /**/ literacy /**/ `tuple`i'' gdppcl, robust cluster(ccode) 
		regsave using /**/ eba_pdest_24 /**/, tstat pval ci autoid append detail(all) 	
						}							 					
																																												
 ***************
 ** 25. lp_catho80 **
 ***************
		
		probit  F.pdest  /**/ lp_catho80  /**/
		regsave using /**/ eba_pdest_25 /**/, replace
		
		tuples 																///
		  agr america_l america_n conflict_ext conflict_int 				///
		 cpi easia_pacific energyuse ethnic_tension ethnic  				///
		 europe execleft execright fdi food      							///
		 homicides infant KOFSoGI language legal_british    				///
		 legal_french legal_german legal_scan literacy      				///
		 /**/  lp_protmg80 media mena migration 							///
		 mortality_under5 lp_muslim80 oecd peacekeepersfl peacekeepersl 				///
		 polity2 pop prevent prostitutionlaw protect 						///
		 pve refugees religion rle sasia 									///
		 ssa tourism trade unemployment urban 								///
		 vae visa wecon  english 										///
		 french german spanish portuguese landlocked 						///	
				, max(3) min(3)  kronecker
		
		forvalues i = 1/`ntuples'   {
		 probit  F.pdest  /**/ lp_catho80 /**/ `tuple`i'' gdppcl, robust cluster(ccode) 
		regsave using /**/ eba_pdest_25 /**/, tstat pval ci autoid append detail(all) 	
						}							 					
																																												
 ***************
 ** 26. lp_protmg80 **
 ***************
		
		probit  F.pdest  /**/ lp_protmg80  /**/
		regsave using /**/ eba_pdest_26 /**/, replace
		
		tuples 																///
		  agr america_l america_n conflict_ext conflict_int 				///
		 cpi easia_pacific energyuse ethnic_tension ethnic  				///
		 europe execleft execright fdi food      							///
		 homicides infant KOFSoGI language legal_british    				///
		 legal_french legal_german legal_scan literacy      				///
		 lp_catho80 /**/   media mena migration 							///
		 mortality_under5 lp_muslim80 oecd peacekeepersfl peacekeepersl 				///
		 polity2 pop prevent prostitutionlaw protect 						///
		 pve refugees religion rle sasia 									///
		 ssa tourism trade unemployment urban 								///
		 vae visa wecon  english 										///
		 french german spanish portuguese landlocked 						///	
				, max(3) min(3)  kronecker
		
		forvalues i = 1/`ntuples'   {
		 probit  F.pdest  /**/ lp_protmg80 /**/ `tuple`i'' gdppcl, robust cluster(ccode) 
		regsave using /**/ eba_pdest_26 /**/, tstat pval ci autoid append detail(all) 	
						}							 					
																						 ***************
 ** 27. media **
 ***************
		
		probit  F.pdest  /**/ media  /**/
		regsave using /**/ eba_pdest_27 /**/, replace
		
		tuples 																///
		  agr america_l america_n conflict_ext conflict_int 				///
		 cpi easia_pacific energyuse ethnic_tension ethnic  				///
		 europe execleft execright fdi food      							///
		 homicides infant KOFSoGI language legal_british    				///
		 legal_french legal_german legal_scan literacy      				///
		 lp_catho80 lp_protmg80 /**/    mena migration 						///
		 mortality_under5 lp_muslim80 oecd peacekeepersfl peacekeepersl 				///
		 polity2 pop prevent prostitutionlaw protect 						///
		 pve refugees religion rle sasia 									///
		 ssa tourism trade unemployment urban 								///
		 vae visa wecon  english 										///
		 french german spanish portuguese landlocked 						///	
				, max(3) min(3)  kronecker
		
		forvalues i = 1/`ntuples'   {
		 probit  F.pdest  /**/ media /**/ `tuple`i'' gdppcl, robust cluster(ccode) 
		regsave using /**/ eba_pdest_27 /**/, tstat pval ci autoid append detail(all) 	
						}							 					
																																												
***************
 ** 28. mena **
 ***************
		
		probit  F.pdest  /**/ mena  /**/
		regsave using /**/ eba_pdest_28 /**/, replace
		
		tuples 																///
		  agr america_l america_n conflict_ext conflict_int 				/// 1-5
		 cpi easia_pacific energyuse ethnic_tension ethnic  				/// 6-10
		 europe execleft execright fdi food      							/// 11-15
		 homicides infant KOFSoGI language legal_british    				/// 16-20
		 legal_french legal_german legal_scan literacy      				/// 21-24
		 lp_catho80 lp_protmg80 media /**/     migration 					/// 25-28
		 mortality_under5 lp_muslim80 oecd peacekeepersfl peacekeepersl 	/// 29-33
		 polity2 pop prevent prostitutionlaw protect 						/// 34-38
		 pve refugees religion rle sasia 									///39-43
		 ssa tourism trade unemployment urban 								///44-48
		 vae visa wecon  english 											///
		 french german spanish portuguese landlocked 						///	
				, max(3) min(3)  kronecker cond(!(20&21) !(20&28) !(21&28))
		
		forvalues i = 1/`ntuples'   {
		 probit  F.pdest  /**/ mena /**/ `tuple`i'' gdppcl, robust cluster(ccode) 
		regsave using /**/ eba_pdest_28 /**/, tstat pval ci autoid append detail(all) 	
						}							 					

***************
 ** 29. migration **
 ***************
		
		probit  F.pdest  /**/ migration  /**/
		regsave using /**/ eba_pdest_29 /**/, replace
		
		tuples 																///
		  agr america_l america_n conflict_ext conflict_int 				/// 1-5
		 cpi easia_pacific energyuse ethnic_tension ethnic  				/// 6-10
		 europe execleft execright fdi food      							/// 11-15
		 homicides infant KOFSoGI language legal_british    				/// 16-20
		 legal_french legal_german legal_scan literacy      				/// 21-24
		 lp_catho80 lp_protmg80 media mena /**/      						/// 25-28
		 mortality_under5 lp_muslim80 oecd peacekeepersfl peacekeepersl 	///
		 polity2 pop prevent prostitutionlaw protect 						///
		 pve refugees religion rle sasia 									///
		 ssa tourism trade unemployment urban 								///
		 vae visa wecon  english 										///
		 french german spanish portuguese landlocked 						///	
				, max(3) min(3)  kronecker cond(!(20&21) !(20&28) !(21&28))
		
		forvalues i = 1/`ntuples'   {
		 probit  F.pdest  /**/ migration /**/ `tuple`i'' gdppcl, robust cluster(ccode) 
		regsave using /**/ eba_pdest_29 /**/, tstat pval ci autoid append detail(all) 	
						}							 					
																																												
***************
 ** 30. mortality_under5 **
 ***************
		
		probit  F.pdest  /**/ mortality_under5  /**/
		regsave using /**/ eba_pdest_30 /**/, replace
		
		tuples 																///
		  agr america_l america_n conflict_ext conflict_int 				///
		 cpi easia_pacific energyuse ethnic_tension ethnic  				///
		 europe execleft execright fdi food      							///
		 homicides infant KOFSoGI language legal_british    				///
		 legal_french legal_german legal_scan literacy      				///
		 lp_catho80 lp_protmg80 media mena migration   						///
		 /**/  lp_muslim80 oecd peacekeepersfl peacekeepersl 					///
		 polity2 pop prevent prostitutionlaw protect 						///
		 pve refugees religion rle sasia 									///
		 ssa tourism trade unemployment urban 								///
		 vae visa wecon  english 										///
		 french german spanish portuguese landlocked 						///	
				, max(3) min(3)  kronecker
		
		forvalues i = 1/`ntuples'   {
		 probit  F.pdest  /**/ mortality_under5 /**/ `tuple`i'' gdppcl, robust cluster(ccode) 
		regsave using /**/ eba_pdest_30 /**/, tstat pval ci autoid append detail(all) 	
						}							 					
																																																			 
						
						
						
						
						
						
	
